The Astrophysical Journal, accepted 

Preprint typeset using I^'T^]X style emulateapj v. 08/29/06 



A MULTIBAND STUDY OF THE GALAXY POPULATIONS OF THE FIRST FOUR 
SUNYAEV-ZELDOVICH EFFECT SELECTED GALAXY CLUSTERS 

A. Zenteno,^'^ J. Song,"' S. Desai,^ "" R. Armstrong,'' J. J. Mohr/'^'' C.-C. Ngeow,' * W. A. Barkhouse; 
S. S. Allam,^ K. Andersson,^ G. Bazin/''-^ B. A. Benson,"' E. Bertin," M. Brodwin, ^^■^'^ 
E. J. Buckley-Geer,^ S. M. Hansen,^-' F. W. High,"' H. Lin,^ Y.-T. Lin,^'-^*^ J. Liu,^'^'" A. Rest,^" 
R. C. Smith, B. Stalder,"^ A. A. Stark, D. L. Tucker,'^ and Y. Yang* 

The Astrophysical Journal, accepted 

ABSTRACT 

We present first results of an examination of the optical properties of the galaxy populations in SZE 
selected galaxy clusters. Using clusters selected by the South Pole Telescope survey and deep multi- 
band optical data from the Blanco Cosmology Survey, we measure the radial profile, the luminosity 
function, the blue fraction and the halo occupation number of the galaxy populations of these four 
clusters with redshifts ranging from 0.3 to 1. Our goal is to understand whether there are differences 
among the galaxy populations of these SZE selected clusters and previously studied clusters selected in 
the optical and the X-ray. The radial distributions of galaxies in the four systems are consistent with 
NEW profiles with a galaxy concentration of 3 to 6. We show that the characteristic luminosities in 
griz bands are consistent with passively evolving populations emerging from a single burst at redshift 
z = 3. The faint end power law slope of the luminosity function is found to be on average a ~ —1.2 
in griz. Halo occupation numbers (to to* + 2) for these systems appear to be consistent with those 
based on X-ray selected clusters. The blue fraction estimated to 0.36L*, for the three lower redshift 
systems, suggests an increase with redshift, although with the current sample the uncertainties are 
still large. Overall, this pilot study of the first four clusters provides no evidence that the galaxy 
populations in these systems differ significantly from those in previously studied cluster populations 
selected in the X-ray or the optical. 

Subject headings: galaxies: clusters - galaxies: evolution - galaxies: formation - cosmology: observa- 
tions 
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1. INTRODUCTION 

Galaxy clusters can be readily discovered or selected 
using optical or IR emission from their member galaxies, 
X-ray emission from the hot intracluster medium and 
now even by the impact of this intracluster medium on 
the cosmic microwave background temperature toward 
these systems. First, from optical observations, |Abell| 
(19581 identified, catalogued and characterized clusters 
of galaxies using classification criteria like compactness^ 
distance, and richness. Later, new optical surveys 
added other optical properties to the clusters. Luminos- 
ity function, radial profile, blue fraction, dwarf-to-giant 
ratio, among others, became tools for understanding dif- 
ferent physical processes in the galaxy cluster environ- 
ment. 

With the advent of space based astronomy new prop- 
erties of clusters of galaxies were discovered. Strong X- 
ray emission made the galaxy clusters some of the most 
luminous objects in the universe, and their properties 
like X-ray luminosity, temperature, and mass have been 
compiled in several X-ray selected cluster surveys (see, 
Giacconi et "aI][T972| |Voges et al.]p99l [2000l [Bohringer 



et al.||2004[ for example) 

In the infra-red regime, the properties of clusters have 
been studied mainly re lying on the X-ray or o ptical clus 
ter identification (see, 



de Propris et al. 



)r or op 
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Lin et al. 



[20031 [20041 [Toft et al.[|20041 [L)e Propris et al.|[W7l 

timore, MD 21218 

21 Cerro Tololo Inter- American Observatory, National Optical 
Astronomy Observatory, La Serena, Chile 



Muzzin et al.||2007a|b| |Roncarelli et al.|[20T0 among 
others) . From IR selected clusters, some of the first stud- 
ies analy zed the cluster populati ons based on individual 
clusters (Stanford et al. 1997, 2005|. Later, systematic 
searches of clusters in the infrared became feasible with 
the operation of space telescopes and with ground based 
telescopes with advanced IR d etectors. S urveys such as 
FLAMEX fElston et al -1120061, UKIDSS ( |van Breukelm 



et al.||2006p , b'L ti jMiIzzin et al. 



Shallow Survey ( Eisenhar^^^^ 



2008t and the IRAC 

have delivered 



cluster catalogs, at high rcdshift, allowing initial system- 
atic characterization of the galaxy populations on those 
systems. 

In the millimeter re gime, the use of the Sunya ev- 
Zel'dovich Eff'ect (SZE, (Sunyaev fc Zerdovich||1972[ ) as 
a selection method for cl uster of galaxies has rec ently 



produced the first r esults ( Staniszewski et al. 2009 Van 



derlinde et al.|2010[ ) . The use of the S/E eft'ect for cluster 
detection has several advantages. A catalog of SZE se- 
lected clusters is approximately mass limited, nearly red- 
shift independent and the o bservable signature is closely 



related to the cluster mass (Birkinshaw 



et al. 2002|, making it less prone to 
selection 



1999 Carlstrom 



e biased in the 



In particular, an SZE selected cluster sam- 
ple provides an opportunity to systematically study the 
galaxy populations and its redshift evolution in clusters 
of the same mass range over a wide range of redshift. 

In this paper we use tools developed for optical stud- 
ies to analyze the galaxy populations of the first four 
SZE selected clusters publi shed by the South Pole T ele- 
scope (SPT) collaboration { Staniszewski et al.|2009 |. As 
well as being among the first SZ selected systems, these 
clusters are among the most well studied. This sample 
has been imaged deeply in th e optical Blanco Cosni ology 
Survey, studied in the X-ray ( Andersson et al.|2010| ), tar 



geted spectroscopically for redshifts ( High et al.||2010[ ) 
and the BCS data have been used to estimate weak lens 

Also these four systems 



ing masses ( Mclnnes et al 
span a broad range in rec 



2009) 

shift and mass, much like the 
larger samples th at have been published s o far ( | Vander 
linde et a l. 20l"0{ I Williamson et aL||2011[ ). In this piloi 
study, we study tne luminosity function, the radial pro- 
file, the Halo Occupation number and the blue fraction, 
in an effort to answer a basic question: Are the galaxy 
populations from these first SZE selected clusters any dif- 
ferent than the populations in clusters selected by other 
means? 

The paper is organized as follows: §2 describes the ob- 
servations and data reduction. In section §3, properties 
of the clusters, such as redshift and mass, are described. 
In §4 we study the galaxy populations in the clusters, 
presenting the main results. Conclusion of this study are 
presented at section §5. Magnitudes are quoted in AB 
system. 

We assume a fiat, ACDM cosmology with Hq = 100ft, 
km s~^ Mpc~^, h — 0.702, and matter density 0„, = 
0.272, according t o WMAP7 
matsu et al.]|2010 |. 
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2. OBSERVATIONS AND DATA REDUCTION 
2.1. Blanco Cosmology Survey 



The Blanco Cosmology Survey^ (BCS) project was 
awarded 60 nights from the NOAO (National Opti- 
cal Astronomy Observatory) survey program starting in 
semester 2005B. Data were gathered in 2005-2008 us- 
ing the Blanco 4-meter telescope located at Cerro Tololo 
Inter-American Observatory^ in Chile. The telescope is 
equipped with a wide field camera called the Mosaic2 im- 
ager, which consists of an array of eight 2K x AK CCDs. 
The pixel scale of Mosaic2 imager is 0.27 arc-second per 
pixel, leading to a field of view of about 0.36 square de- 
gree. The observations were carried out to obtain a deep, 
four band photometric survey {g, r, i and z) of two 50 
deg^ patches of the southern sky centered at 23h00m,- 
55° 12" and 05h30m,-55°47". These regions were chosen 
to enable observations by three mm-wavelength survey 
experiments (the SPT, the Atacama Pathfinder Exper- 
iment (APEX) and the Atacama Cosmology Telescope 
(ACT) experiments) . On photometric nights we also ob- 
served several standard star fields that contain stars with 
known magnitudes. This approach allows the calibration 
of our data to the standard magnitude system. In ad- 
dition, we obtained deep imaging of several fields over- 
lapping published spectroscopic surveys to enable cali- 
bration of photo-z's using samples of many thousands of 
spectroscopic redshifts. The data volume we collected for 
the BCS observation was about 20 to 30 Gigabytes/night. 

The first three seasons (2005 to 2007) of the BCS imag- 
ing data were processed in 2008 and 2009 using version 
3 of the data management system developed for the up- 
coming Dark Energy Survey (DES). Details of the DES 
data manag ement system can b e found in [Ngeow et al.| 
((20061) and |Mohr et alT] ( |2008[ ). A brief description is 
presented in this paper! Data parallel processing was 
carried out primarily on NCSA's TeraGrid IA-64 Linux 
cluster. The pipeline processing middleware developed 
within the DES data management system provides the 
infrastructure for the automated and robust execution of 
our parallel pipeline processing on the TeraGrid cluster. 

To remove the instrumental signatures, the raw BCS 
images were processed using the following corrections: 
crosstalk correction, overscan correction, bias subtrac- 
tion, flat fielding, fringe and illumination correction. Bad 
columns and pixels, saturated pixels and bright star ha- 
los, and bleed trails are masked automatically. Wide 
field imagers have field distortions that generally deviate 
significantly from a simple tangent plane, and there are 
typically telescope pointing erro rs as well. The AstrO- 
matic code SCAMP ( | Bert in 20061 was used to refine the 
astrometric solution by matching the detected stars in 
BCS images to the USNO-B catalog. We adopted the 
PV distortion model that maps detector coordinates to 
sky coordinate using a third order polynomial expansion 
of distortions, across each CCD, relative to a tangent 
plane. The DES data management system is using an 
experimental version of t he AstrOmatic tool SExtractor 
( [Bertin fc Arnouts||1996| ). This experimental version in- 
cludes model fitting photometry and improved modes of 
star-galaxy classification to detect and catalog astronom- 



http : //cosmology . Illinois . edu/BCS/ 

Cerro Tololo Inter-American Observatory (CTIO) 
is a division of the U.S. National Optical Astronomy 
Observatory (NOAO) , which is operated by the Association 
of Universities for Research in Astronomy (AURA) , under 
contract with the National Science Foundation. 
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ical objects in the images. We harvested a wide range of 
photometric and astrometric measurements (and their 
uncertainties) for each object during this cataloging. 

For the photometric nights that include observation of 
the standard star fields, we determined the band depen- 
dent (atmospheric) extinction coefficients (fc) together 
with CCD and band dependent photometric zeropoints 
(a) and instrumental color terms (6). Specifically, the 
equation we constructed for each star in the standard star 
fields is minst-rristd = Si ^ [ai + hi{l:\C)]+kX , where 
Wj = 1 if the standard star is on CCD i\ Wi = other- 
wise. In this equation, niinst and mstd are the instru- 
mental and the true magnitudes for the standard stars, 
respectively, AC is the color offset of the standard stars 
from a reference color, and X is the airmass. The stan- 
dard star fields include the SDSS Stripe 82 fields and the 
Southern Standard Stars Network fields^. The resulting 
photometric solutions were then used to calibrate the 
magnitudes for other astronomical objects observed on 
the same night. 

The nightly reduced and astrometric refined images 
were remapped and coadded to a pre-defined grid of tiles 
(which is a rectangular tangent plane projection, with 
~ 36 arc-minute on a side (hereafter the BCS t iles) in 
the sk y using another AstrOmatic tool, SWarp (Bertin 
et al.|[2002) . During this coaddition we carry out a PSF 
homogenization across each tile and within each band to 
match the PSF to median delivered seeing in that part 
of the sky. The zeropoints for the flux scales for these in- 
put remap images are determined using different sources 
of photometric information, including direct photomet- 
ric zeropoints which are derived from the photometric 
solution on photometric nights, relative photometric ze- 
ropoints determined using all pairs of images that over- 
lap on the sky and the color behavior of the stellar locus 



( High et ah] 2009 1 . We determine the zeropoints for all 
images by doing a least squares solution using the con- 
straints described above. During co-addition, we use a 
weighted mean combine option in SWARP. The coadded 
images are b uilt in each band fo r a given coadd tile, then 
a image (Szalay et al. 1999) is created for detection 
and cataloging to ensure each object will have measure- 
ments in the griz bands extracted from the same portion 
of the object. 

2.2. Completeness 

For this work we estimate the completeness of the 
BCS tiles from the comparison of their griz source 
count histograms and those extracted from the deeper 
Canada- Fr ance-Hawaii- Telescope L egacy Survey survey 



(CFHTLS, jBrimioulle et~aL][2008| private communica- 



tion) . Specifically, we used count histograms from the D- 
1 1 sqr. degree patch at high galactic latitude (1= 172.0°; 
b = -58.0°) from the CFHTLS Deep Field, whose mag- 
nitude limit is beyond r=27 and the seeing is better than 
1.0" and 0.9" for g and r-zz^respectively Dividing both 
count histograms (see Fig.nJ we can estimate the level of 
completeness in the different tiles in each band. We can 
use this completeness estimate for each field to account 
for the missing objects as we approach the full depth of 
the photometry. Table[l]contains the magnitude limits in 
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Fig. 1. — We estimate the completeness of our BCS coadds 
by comparing objects counts from them with counts from deeper 
CFHT data. The resuhing completeness curve is fitted by an error 
function, which later is used to correct for the missing galaxies and 
to define 90% and 50% completeness limits for analysis. Here is an 
example for the SPT-CL J0516-5430 field. 



each band corresponding to 50% and 90% completeness 
for the tiles used in our analysis. 

3. BASIC PROPERTIES OF THESE SPT CLUSTERS 

The basic properties of these SPT selected clusters, 
including the characteristics of the o ptical counterparts 
are presented in Staniszewski et al. (200 9J and are fur 



ther discussed in follow on papers (Mcnant eau et al.|2009 
Mclnnes et"aLl[2009l [High et al.[ 2010; A ndersson et al". 
20101. Several spectroscopic redshifts are now available 



as well as Chandra X-ray observations, providing dra- 
matically improved mass information which enables the 
kind of galaxy population study we undertake here. De- 
spite being a small sample, these four clusters are among 
the most well studied SZ selected systems and their mass 
and redshift distributions range are similar to the whole 
SPT published cluster sample (see Fig.|2|. In particular, 
these masses and redshifts are used to estimate the pro- 
jected cluster virial radius in which the optical properties 
are measured. These properties are presented below. 

3.1. Redshifts 

The spectroscopic redshifts of the four cluster are now 
available (Table[2|. SPT-CL J0516-5430, a known cluster 
identified in the Ab ell supplementary southern catalog 
(AS0520, lAbell et al. 1989) and in the REFLEX survey 
(RXC J0516.65430, Bohringer et al. 2004), had a redshift 
of 0.2 94 (Guzzo et al.,1999) and 0.2952 ( .Bohringer et al. 
2004[ ), values obtained using 8 galaxy spectra from the 
ESU Key Programme. 

For the other three clusters, spectroscopic data has re- 
cently be en acquired. Usin g LDSS-3 on Magellan Clay 
telescope. High et al. (20101 reported redshifts of 0.7648 
for SPT-CL 0528-530ini5rd 0.4 626 for SPT-CL 0509- 
5342. For SPT-CL J0546-5345, iBrodwin et all ([20101 



used IMACS on Baade Magellan telescope to measure a 
redshift of 1.0665. 



http : / /wm-stai . f nal . KOv/Southern_ugriz/index . html 




Details can be found at 


http : / /www . ast .obs-mip.fr/ 


article212.html 



4 



TABLE 1 

Completeness limits for each tile for each filter for 90%/50% completeness 



ID 


R.A. 


decl. 


g 


r 


i 


z 




[deg] 


[deg] 


500 sec 


600 sec 


1350 sec 


705 sec 


SPT-CL J0516-5430 


79.15569 


-54.50062 


23.18/24.24 


22.73/23.87 


22.20/23.47 


21.87/23.19 


SPT-CL J0509-5342 


77.33908 


-53.70351 


23.72/24.78 


23.29/24.51 


23.10/24.23 


22.45/23.78 


SPT-CL J0528-5300 


82.02212 


-52.99818 


23.70/24.62 


23.42/24.32 


22.93/23.94 


22.23/23.37 


SPT-CL J0546-5345 


86.65700 


-53.75861 


23.34/24.31 


22.87/23.90 


22.48/23.64 


21.97/23.08 



Wll sample ^ 
VIO+AIO sample ; 
SQ9 sample ; 



0.2 0.4 0.6 O.S 1 1.2 

redshift 

Fig. 2. — Masses and redshifts of th e SPT cluster samp l e pub - 
lished to date. Open c ircles are from 'Williamson et al.' ('201 1|, 
open triang les are from [And crsson ot al. (2010) and Vandcrlinae| 
|et al.| ( [2010^ , and filled squares trom ^ijtaniszcwski ct al._ ( .2U0 9>, the 
sample here studied. The M2oo,crit mass estimations come from 
X-ray observations where that is possible, or from the SPT detec- 
tion significance. In the latter case, masses have been converted 
from their native estimation M200,mean to M2oo,crit assuming a 
concentration parameter of c = 5 for the halo mass (see tablel2]l 
under assumed cosmology. It can be seen that our sample spans 
on the redshift and mass space for the latest SPT sample. 



3.2. Cluster masses 

As mentioned in S|3j the optical analyses performed in 
this work require an estimate of the projected cluster 
virial radius. For this purpose, along with spectroscopic 
redshifts, X-ray masses estimations are used. We adopt 
mass estimates defined with respect to the critical den- 
sity. 

As it has been previously mentioned, SPT-CL 0516- 
5430 is a previously known cluster, and its mass has also 
been estimated. With the name of RXC J0516. 6-5430 



in the REFLEX survey, [Zhang et al.| ( |2006[ ) used XMM- 
Newton to find a M500 of (674 ± 2.1)IlT'* Mq. Also, 
recent X-ray observations of the four clusters have been 
performed, and the mass estimation of SPT-CL 0517- 

5430 has been refined. 

Usin g Chandra and XMM-Newton, Andersson et al.] 
(2010) reported X-ray measur ements of 15 of the 21 SZE 
selected clusters presented in Vanderlinde et al. (2010). 
The observations of those clusters, which include the 
original first four clusters, have been designed to deliver 
around 1500 photons within 0.5r5oo, in order to enable 
measurement of the ICM mass and ICM temperature, 
allowing a mass estimation t hrough a M500 — Yx scaling 
relation ( Vikhlinin et al.||2009j ) with approximately 15% 
accuracy. 



From X-ray Mspn .yy and spectroscopic redshifts An- 
dersson et al. (2010) estimated the physical r^ooXx' both 
dehned with respect to the critical density. Here, those 
r500,Yx/M5oo,Yx are transformed tor^ooji/MafloJi- us- 



ing the Navarro, Frenk, & White ( Navarro et al.||1997 
hereafter 'NEW' profile) radial mass prohle with con 
centration of 5 for the dark matter halo, which implies 

r200,Yx = 1-51 X ^500, / -^200, Yx = 1-38 X MsqcYx Con- 
version. The angular projection is calculated using the 

spectroscopic redshifts. M2oo,Yx "^^U as r2oo,Yx ^-re 
listed in Table [2 

4. CLUSTER GALAXY POPULATIONS 

The galaxy populations in clusters have been studied 
using several techniques and selection processes. Clusters 
of ga laxies have been selected mainly from optical im- 



2007 



ages ( |Abell||1958[ [Abell et aL||1989| [Koester et al 
for exam ple) an d through their X-ray emission ( Ebel 
ing et ar]|1996 Vikhlinin et al. 1998 |B6hringer et al 
2004[ aniong others). A selection of clusters based on 



their Sunyaev-Zeldovich signature promises a less biased 
selection method, as it is likely to be less affec ted by 
projections or false clusters like optical surveys ( Lucey 
T9831 |Sutherland|[T988l |Collins et ai:][T995l |Cohn et al. 
2007|), and its mass selection function is nearly redshift 
independent, making the cluster sample more homoge- 
neous than X-ray surveys in redshift space. Also, the 
SPT survey will be able to find the most massive clus- 



ters o f galaxies in the universe (Carlstrom et al. 2002 



2009[ ). Once completed, the size, redshift extent and the 
degree of completeness of the SZE selected cluster sam- 
ple will be ideal for statistical analysis of astrophysical 
properties in high mass clusters. Here we focus on the 
first four SPT selected clusters, which are all high mass 
systems extending over a broad redshift range. 

4.1. Radial distribution of galaxies 

The radial distribution of galaxies in clusters can be 
used to further our understanding of the cluster envi- 
ronment physics. For example, from N-body and gas 
dynamical simulations, which include radiative cooling, 
star formation, SN feedback, UV heating, etc., |Nagai fc 



Kravtsov ( 2005 ) produced radial distributions consistent 
with observations of X -ray selected cluster s a mples from 
Carlberg et" aLl (|1997[) and [Lin et aI 1(|2004[). |Saro et al 
( |2006p , using hydrodynamical simulations, also showed 
an agreement between the radial distribution of the sim- 
ulate d galaxies and X-ray a nd optically selected clusters 



from Popesso et al.| ( |2007a ) 

For the following analysis we define the cluster cen- 
ter to be the position of the observed brightest clus- 
ter galaxy member (BCG; coordinates are listed in Ta- 



blejlj), which agree with its X-ray center (Andersson et al 
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TABLE 2 

X-RAY MASSES, SPECTROSCOPIC REDSHIFTS AND CLUSTER PARAMETERS. 



ID 



[Mpc] 



T20a,Yx 

[arcminl 



S/N 



^red gal 



^all gal 



8.34 
4.54 
2.61 
2.57 



9.42 
6.61 
5.45 
7.69 



fO.81 
-0.73 
3^^g-f3.50 



4.65: 



Note. 



SPT-CL J0516-5430 16.38 ± 1.72 0.2952^ 2.21 

SPT-CL J0509-5342 7.51 ± 0.83 0.4626^ 1.61 

SPT-CL J0528-5300 4.11 ± 1.19 0.7648" 1.17 

SPT-CL J0546-5345 7.37 ±0.85 1.0665'' 1.27 

"1.38 X M500 from Andersson et al. (20101, assuming a concentration parameter of 



2.79: 



1.39 

1.98 
1.37 



4.02^ 



'0.52 

1 04+7.44 

3.23+^-^'^ 
4.04 



0.55 

1.92 
1.31 



(201(1 



1.51 X r^oo^Yx from Andersson et al ^ . 

measured in 150 GHz SPT maps from from [Vanderlinde et aL (20101 
red ga laxies 
(20041) 3 



5 for the mass halo. 

assuming a concentr ation parameter of c = 5 for the mass halo. "^The S/N 



Concentration parameter from the N FW fitting of the 

Concentration param eter from the NFW fitting of the all galaxies. -^S pectroscopic redshift fr om Bohringer et al. 

Spectroscopic redshift from High et al.| ( |2010| . ''Spectroscopic redshift from Brodwin et al.| (2010l. 




SPT-CL 0509-5342 



SPT-CL 0528-5300 
cg= 5.93 .2 jj 




0.4 0.8 1.2 1.6 2 2.4 2.8 0.4 0.8 1.2 1.6 2 2.4 2.8 
^200 ^200 

Fig. 3. — Radial profiles for the red galaxy population to 0.36L* 
(same as blue fraction), binned to obtain similar signal to noise. 
These profiles are centered on the BCG and extend to 3r2oo to allow 
the background and the cluster profile to be fit simultaneously. All 
radial profiles are consistent with NFW profiles with concentrations 
c ~ 4. 



2010 1. In order to compare different studies we estimate 
the concentration parameter from the NFW surface den- 
sity profile. We obtain the NFW surface density by in- 
tegrating the three-dimensional number density profile 
n(x) = nQX ~^(l+x)^^ along the line-of-sight (seelBartel- 



mann|1996 1 where x = Cgr/r2oo, is the normalization 
and Cg is the concentration parameter. On stacked clus- 
ter data it is customary to fit both parameters, no and Cg. 
In our case the NFW fit is done over single cluster data 
to a common magnitude limit, and this leads to consider- 
able uncertainties in the parameters of the NFW profile. 
In order to minimize this problem, we introduce the ob- 
served number of galaxies in the equation. Integrating 
the NFW surface density over the projected area we can 
derive uq — no{nobs,Cg) allowing us to fit the NFW den- 
sity profile as a single parameter function. Also, after a 
statistical background correction, a background fitting is 
performed along with the NFW fit. Such background fit 
is limited within the Poisson uncertainty of the observed 
background. 

The radial surface density profiles are constructed us- 
ing both the red±blue and the red galaxy population 
defined from the col or-magnitude diagram of the red se- 
quence (see sec. 4.4). The galaxy population is also se- 
lected performing a cut in brightness, selecting galaxies 



which are fainter than the BCG and brighter than a com- 
mon limit of 0.36L*. Th e error bars a re computed using 
small number statistics ( Gehrels|1986 ). The background 
is statistically subtracted and a second correction is ap- 
plied fitting it to a radius of 3r2oo- Finally the data is 
presented using radial bins of constant signal-to-noise of 
3.5 (see Fig. [3]). 

Some corrections are applied to these profiles. In the 
area calculation for each radial bin in Fig. [3j the area cov- 
ered by saturated stars was excluded in order to avoid 
an under estimation of the surface density. This is espe- 
cially important in the case of SPT-CL J0509-5342 where 
several bright stars close to the BCG are blocking the de- 
tection of galaxy cluster members, covering about 50% 
of the area at a radius of 0.4r2oo- 

The concentration parameters found are shown in Ta- 
ble [2] With a concentration in the range of c « 3 — 6, 
the clusters agree at la confidence. We note that the 
blue+red distribution tends to be less concentrated than 
the red population alone, which is consistent with pre- 
vious analyses where a higher concentratio n is seen in 



the red population (e.g. Goto et al. 2004). The con 



centration we find is in agreement with concentrations 
drawn from X-rav selected c lusters of galaxies. For ex- 
ample, [Carlberg et"al| \l997\ found a Cg of 3.70^?;^^ at 
95% confidence, using 16 clusters from the CNOC survey 
with a median redshfit of ~ 0.3 for a similar n iass range 



(2 X X AfiS; Carlberg et al](1996)). 



(20041, from stacked 2MASS K-band data on 93 nearby 



Lin et al. 



X-ray selected clusters, found a value of Cg — 2.90^022 
in a wider 3 x Mq-2 x Mq range. Both are consistent 
with our results. 

We also found agreement between our concentration 
parameter and the conce ntration parameter found for 



optical selected clusters. Biviano & Poggianti (20091 



found, studying 19 intermediate redshift (0.4 < z < 0.8 
0.7 < M200 < 13.6 X lO^^Mf^) EDisCS±MORPHS clus 



'200 rj-u ^ ru J^-'o 
ters, a concentration parameter of c 



3.2 



+4.6 



Also, 



Johnston et al. (20071 found, using the SDSS sample, a 
concentration parameter of C2oo|i4 = 4.1±0.2siaf ±1.23j,s 
for a cluster mass of M = 10^'^h~^ Mq. In summary, we 
find no evidence that SZE selected clusters exhibit differ- 
ent galaxy radial distributions than in optical and X-ray 
selected clusters. 
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Fig. 4. — Here we plot the (with 1 a uncertainties) for each 
band that results from Schechter function fits to the luminosity 
function with free parameters m*, (p* and ct (a fixed where noted). 
We limit the range of m* to be fainter than the identified BCG 
for each cluster. The continuous lin e represents a passiv e evolu tion 
single burst model at 2 = 3 from pBruzual &; Chariot] | |2003] |. It 
is clear that these SPT selected clusters have galaxy populations 
consistent with this simple evolution model. 



4.2. Luminosity functions 

The luminosity function (LF) is an important tool for 
testing theories of galaxy formation and evolution. For 
example, ever more complex simulations can be tested 
against the LF, as an observational constraint, to probe 
our understanding of th e evolutio n of gala xies in the clus - 
ter's environment (Romeo et al.||2QQ5; Saro "eral|[2006| . 
With clusters of similar masses we can study the Li"' as a 
function of redshift and with the LF parameters we can 
calculate the Halo Occup ation Number (HON) and test 



the N-M scaling relation (Lin et al 
The LF can be descr ibed by t 



2004 2006) 



Schechter function ( Schechter||1976 ), 



le three parameter 



0(m) = 0.41n(10) $*10"-'*'"*-™'("+i) exp(-10"-'*('"*-'")) 

where $* is the normalization, m* is the characteristic 
magnitude and a accounts for the faint end power law 
behavior of the function. 

The construction of the LF is done assuming that the 
observed LF in the cluster area is the superposition of the 
cluster LF and the background/foreground non-cluster 
LF. To recover the cluster LF we subtract the galaxy 
source count, rescaled by the area, from the observed 
LF. Given the wide range in redshift we present the LF 
in the four griz bands. 

The are a of the cluster is defined by our estimation of 
^"200 (see ^3.2 and Table [2|, and the area of the back- 
ground is the tile area (3(^x36') minus the cluster area. 
The bright end limit of the LF is defined by the clus- 
ter's BCG while the faint end limit is defined by its 
completeness at 90% or 50%, depending on the redshift, 
in each band. Below the 100% completeness, the 0.5 
mag bins are corrected using the error function fitted to 
the BCS / CFHT comparison galaxy count histograms de- 
scribed in §2.2| and shown in Fig. ij Finally, the number 
of galaxies, background correctecff is divided by cluster 
volume (in Mpc"^) and the uncertainty is assumed Pois- 
sonian in the total number of galaxies (cluster plus back- 



ground) . 

Below we extract m* and a from our cluster sample 
and compare them to previous results drawn from X-ray 
and optical selected clusters of galaxies. 

4.2.1. Evolution of m* 

Studies of m* evolution in clusters have been done in 
different wavelengths and with different selection meth- 
ods. These studies indicate that the stellar populations 
in many of the cluster galaxies have e volved passively 



after forming at high redsh i ft (see, e.g., Gladders et al 

[19981 |De Lucia et al.||2004l [Holden et aL, 2004 

ct al. 2008j and references therein) . There are several 
indications that m* evolution can be described from a 
single stellar populations (SSP) synthesis model as op- 
tical and X-ray selected clusters. All four of these SPT 
selected clusters in our study show red sequences (see 
Fig. 11), and their color evolution is consistent with col- 
ors derived from a single stellar population (SSP) syn- 
thesis model. 

In order to perform a direct comparison of the bright- 
ness and evolution of the characteristic magnitude, we 
let all the LF variables vary, where possible, and com- 
pare TO* in griz bands derived from the LF fitting to 
that based on the SSP model. The SSP model we use for 
the red galaxy population is constru cted using a Bruzual 
fc Ch ariot synthesis model (BC03; Bruzual & Chariot 



20031 for the red galaxy populations, assummg a sm- 
gle burst of star formation at z=3 followed by passive 
evolution to z=0. We use six different models with six 
distinct metallicities to match the tilt of the color mag- 
nitude relation at low redshift, and we add scatter in the 
metallicity-luminosity relation to reproduce the intrinsic 
scatter in the color-magnitude relation. These models are 
then calibrated, using 51 X-ray clusters that have avail- 
able SDSS magnitudes drawn from the DR7 database. 
Details of the model used can be found in Song et al 
(submitted). As shown in Fig. [4j the SSP model and to* 
in each band are in good agreement, showing that the 
SSP model is an appropriate description of both the col- 
ors and the magnitudes of the more evolved early type 
galaxies in this sample of SZE selected clusters. We will 
use this agreement to carry out a more constrained study 
of the luminosity function. 

4.2.2. Faint end slope 

To learn about the a behavior we take advantage of the 
agreement shown in § [4.2.1 [ between SSP model and the 
data. We adopt m* from the model (see Table [3]) and fit 
for $* and a for each cluster individually. The study of 
the faint end slope a provides us with information about 
the faint galaxy populations in the cluster with respect to 
the more evolved bright end, which is dominated by lu- 
minous early type galaxies. This relation gives us insight 
into competing processes in the hierarchical structure for- 
mation scenario, including the accretion of faint galax- 
ies by the cluster, causing a steep a, and the evolution 
of galaxies inside the cluster through galaxy merging, 
dynamical friction, star formation quenching and other 
processes. 



Using the amoeba simplex minimization routine ( Press 



et al. 1992), <I>* and a are chi-square fitted, and their 



uncertainties are determined by gridding in parameter 



7 
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It = 4.35 Mpc " 
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* =3.45 Mpc' 
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Fig. 5. — Here we plot the luminosity function with best fit 
Schechter function for each band in SPT-CL J0516-5430. Note 
that the BCG had been removed. Best fit parameters are shown 
on the figure, while Table |3] includes best fit and Icr uncertainties. 
100% and 90% completeness limits are noted with vertical dotted 
lines in each panel. 



iUli 



and their contour 



space (see the LF in Fig.^ 
confidence regions at Fig. 

From the hterature, we nnd that our average a ~ - 
1.2 is in agreement at the Icr level with previous studies 
which used samples constructed with different selection 
methods. Fo r example, from an optical work, De Propris 
et al.| poos]) used 60 clusters at z < 0.11 from 2dFGRS 
in the bj band finding a = —1.28 ± 0.03. [Paolillo et al. 



(|2001|) found, on a composite LF of 39 Abell clusters, an 
+"■09 and 



a of -1.07 



+0.09 
-0.07' 



-1.11" 



0.07 



-1.09tEJ ii for Gunn g, 



r and i respectively. 

From X-ray selected samples |Lin et al.| ( |2004[ ) created 
a composite K-band LF of 93 clusters, finding that the 
faint-end slope is well fitted by —1.1 < a < —0.84 in 
agreement w ith our findings within the errors. [Popesso] 
et al. ( 2005 1 , using 97 X-ray selected clusters with SDSS 
photometry, for a redshift z < 0.25, found that a better 
representation of the data is given by two Schechter func- 
tions, characterized by a bright and a faint end slope. 
Comparing to the bright end of the double Schechter 
function with local background subtraction (which is the 
most similar case), the bright end slope, in 1 Mpc ft,^^, 
has aslope a of -1.23±0.11, -1.05±0.13, -1.17±0.13, 
and — 1.06 ±0.12 in r, i, and z, respectively, also agree- 
ing with our findings at th e la level. 



From IR selected clusters Muzzin et al. ( 2008 1 detected 
99 clusters and groups of galaxies and constructed the 
LF in 3.6/im, 4.5/im, 5.8/im, and 8.0/im. Although the 
3.6/xm band is redward of our griz photometry, the LF 
constructed seems to be consistent with a ^ —1. 

The agreement found between the multiband LF pa- 
rameters calculated for our SZE selected clusters, and 
previous studies of galaxy cluster LFs indicates that the 
galaxy populations in these SZE selected clusters are not 
very different from those in clusters selected by other 
means. 

4.3. Halo Occupation Number 



Based on the Press & Schechter formalism (Press &: 
Schechter|1974 ) , the halo occupation distribution (HOU) 



Fig. 6.— Same as Figure[5]but for SPT-CL J0509-5342. 
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Fig. 7.— Same as Figurejslbut for SPT-CL J0528-5300 with the 
100% and 50% completenesslimits noted with vertical dotted lines 
in each panel. 
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FlG. 8.— Same as Figure [7] for SPT-CL J0546-5345 with only i 
and z bands. 



is a powerful analytical tool for understand ing the phys 



ical processes drivin g 



Berlind et al. 2003) 



galaxy formation 
Also the HOD 



(Seljak 



_ 2000 
uses' 



to 



constram cosmological parameters ( [Zheng &: Weinberg 
20071. 
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TABLE 3 

HON PARAMETERS 



ID 




^9 


m; 


HON 

[m; + 2) 






Q 




m* 


HON 

{m* + 2) 




SPT-CL J0516-5430 


3.95+^,11 


-1 29+'^-^'* 


20.73 


327+i^^ 


0.44 


q 70 + 0. 81 
0. / Q 'JQ 


-1.27 


+0.10 
^0 


19.28 


308+|"3^ 


0.81 


SPT-CL J0509-5342 




-1.1" 


22.27 


19111° 


0.35 


4.36l}-^^ 


-1.11 


+0'26 
-0.28 


20.74 




0.58 


SPT-CL J0528-5300 






24.34 






6.8911-JI 


-1.82 


+0.66 
-0.91 


22.52 


149t«^^ 


0.16 


SPT-CL J0546-5345 






25.92 












23.85 






ID 


$* 




m* 


HON 


x' 




Q 


z 


m* 


HON 


x' 










(m* 2) 












{ml + 2) 




SPT-CL JOS 16-5430 


o =0+0.78 
o.oy_Q 74 


1 17+0.10 

J^-J^'-o.io 


18.76 


266+?^ 


1.67 


n (-9 + 0. 82 
'^•'-'^-0.77 


-1.11 


+0.11 
-0.11 


18.42 


248+^5* 


0.93 


SPT-CL J0509-5342 




1 10+0. 27 
^•^■^-0.28 


19.99 




0.39 


3-45ll:f8 


-1.23 


+0.24 
-0.26 


19.57 




0.29 


SPT-CL J0528-5300 


7 lt:+4.57 


1 90+0.46 
J^-^'5-0.52 


21.42 


84tjf 


0.64 


^•0*_3.46 


-1.47 


+0.50 
-0.82 


20.78 


cc;+265 


0.51 


SPT-CL J0546-5345 


25.43ll;°i 


-i.r 


22.83 


337if3 


0.88 


1 r 9+5.57 


-1.1" 


21.75 


2021^1 


1.19 



Note. — "q set to fixed value, m* from model of passive evolution. 
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Fig. 9. — We plot the 68% confidence region for the LF param- 
eters for each cluster and band combination. Panels are arranged 
by band with confidence regions for each cluster where a fit for a 
and <I>* was possible. The current data suggest steeper than nor- 
mal faint end parameters a in two of the clusters and there is a 
tendency for the higher redshift systems to have higher character- 
istic galaxy densities, as expected in an evolving universe. Values 
including those at fixed a are given in Table [3] 



One of the key ingredients in the HOD formahsm is 
{N){M), the mean number of galaxies per halo or the 
Halo Occupation Number (HON) . In the hierarchical sce- 
nario, the HON is expected to increase slower than the 
mass. While a fraction of the accreted galaxies merged, 
the galaxy production becomes less efficient as larger 
haloes are also h otter and less efficient in gas cooling 
( Cole et al.||2000 ). Observationally, several studies with 
cluster samples selected optically and through their X- 
ray emission have been performed, reinforcing that pic- 
ture. For example, from samples of opti c ally s elected 
clusters and groups, Marinoni & Hudson (20021 found 
N oc M°-S3±o.i5 f systems with M > 10^, 
Also, Muzzin et al. 
in the ~ 2 X WMp 



systems with 
(|2007cj) found 
- 2 X IQi^MfT, 



M > IO^^Vs^Mq. 

N /i#0.71±0.11 
-'^'500 

mass range . In the 



X-ray selection method counterpart, Lin et al. (20041 
found, from a sample of 93 nearby clusters and groups, 
N IX. jVf o .87±o.04^ Combining X- rav and optically selected 



clusters, 



Popesso et al. (|2007a 
ound 



A similar picture was 



found N (X Ml 



0.92±0.03 
2fiO, 



by Rines et al. (2004) 



who used nearby X-ray luminous Abell clusters of mass 
- 3 X lO"/i-iM0 and found cx MO-^^io.is^ 

Here we test whether the HON of SZE selected clusters 
exhibits a N (x M^, with (3 < 1, behavior shown by other 
selection methods. 

Due to the small sample presented here, our approach 
is to construct the HON and compare our results to 
the N-M scaling rela t ion and evo lution constraints ob- 



tained by |Lin et al.| (|2004[ |2006[). That scaling rela- 



tion is appropriate in this analysis as it covers the mass 
and redshift range of this SZE sample. The scaling re- 
lation was constructed using X-ray selected clusters in 
the 3 X IO^^Mq - 2 X IO^^Mq mass range using ne arby 
cluste rs with 2MASS K-band photometry, and later, |Lin| 
et al.| (^006), counting galaxies to the depth m*+2, ex- 
panded the study to the 0-0.9 redshift range showing that 
the rel ation does not s trongly evolve. 

( |2004[ ) N-M relation is, 



The Lin et al. 



(36±3)(M2oo/(lO^*/ifoM0)) 



0.87±0.04 



To calculate A'200 we integrate the cluster luminosity 
function to L{m*j^Qjj^j^ + 2) using the parameters of the 
Schechter luminosity function fit, 0*, L* and a computed 
in §4.2[ The total number of galaxies is 



N + with = V(f>^ 



y e 



dy 



where the 1 comes from the BCG, which is not part of the 



LF fitting, V is the cluster volume, and yio 



We use the derived M- 



2 00. 



masses and uncertainties as 



explained in ^3.2 from Andersson et al. (2010) Chandra 
and XMM observations. I'he uncertainty in A200 is es- 
timated by propagating the 1 a uncertainty in and a 
through the integration of the LF to m* + 2. 

The A200 with their X-ray mass for the four clusters 
in the fo ur observed band , along with the HON relation 
found by Lin et al. (2004), are shown in Fig. 



10 Agree- 



ment between these SFT clusters and the published re- 
sults on the X-ray selected sample is good. As with the 
concentration and the LF faint end, there is no signifi- 
cant evidence that the galaxy properties differ from those 
already extracted from previous X-ray selected cluster 
samples. 

4.4. Blue fractions 
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Fig. 10. — We present the halo occupation number 
(HON(m*+2)) within each band for each cluster where the LF 
is measured in more than three bins. Masses and uncertainties on 
the horizontal axis come from X-ray analysis of Chandra observa- 
tions l^ndcrsson et al. 2010). HON uncertainties are derived from 
the variation of HON due to the la uncertainty in the LF (a and 
<&*). The dotted lines show the HON derived fro m a K-band anal- 
ysis o f a large sample of X-ray selected clusters | |Lin et al.||2004] 
|2006| . These SPT selected clusters appear to be neither richer nor 
poorer. 



Another property of the galaxy populations used to 
study their evolution in clusters of g alaxies is the blue 
fraction (f, 



ler 




1984 


tical 


iy se 



fb). In their seminal work Butcher & Oem- 
BO hereafter), using a samples of 33 op- 
ected clusters of galaxies, estimated ft and 
showed that it increased with look-back time (termed 
the Butcher-Oemler e ffect). Later studies , such aslRakos 
k Sch ombert' (1995) (0 < z < 1) and Mar goniner &; 
arCar valho (2000) (0.03< z <0.38), using optically se- 
lected clusters, also have found a strong increase in /t, 
with redshift. 

With the advent of new optical surveys with hundreds 
or thousands of clusters the analyses have been strength- 
ened statistically. Using a sample of salOOO clusters, in 
a wide redshift range (0< z <0. 9), drawn fr om the Red- 



Sequence Cluster Survey (RCS), |Loh et al.|p008J found 



a mild correlation between the red fraction and redshift 



Hansen et al. 



( 2009| ), using thousands of clusters and 

groups from SDSS, found an evolving ft in the two red- 
shift bins studied (0.1-0.25 and 0.25-0.3), also noticing 
that fb evolution was weaker for optical masses above 
10"/i-iMo. 

Studies using samp les of X-ray selected cluste rs, have 
been contradictory. Kodama & Bower (2001) used a 
sample of seven clusters, in the redshift range of 0.23- 
0.43, and found a blue fraction trend consistent with BO, 



while Fairley et al. (2002), using a sample of eight clus- 



ters in a 0.23-0.58 redshift ra nge found virtually no trend 
with redshift. More recently, " 
CFHT MegaCam g and r p 



Urquhart et al.| ( 2010[ ) used 
lotometry on 34 X-ray se- 
lected clusters in the redshift range 0.15-0.41 to study fb 
correlation with other intrinsic cluster properties, found 
that fb correlated with mass {Tx) and redshift. 

Also there are en vironmental factors to be considered. 
Small et al. ( |1998 ) used 10 X-ray selected clusters at 
similar redshiit (0.22 to 0.28) and found a low blue frac- 
tion of 7b = 0.04 ± 0.02 with a variation of Afb = 0.06, 
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Fig. 11. — Color magnitude diagram for galaxies around each 
cluster. The blue population is defined to be more than 0.22 mags 
bluer than the red sequence. Selection in magnitude uses the BCG 
on the bright end and m* -f (5 on the faint end, where this limit 
corresponds to the 90% completeness limit for SPT-CL J0528-5300 
(0.36-LJ^Qjjgj^). The visually identified red sequence cluster galax- 
ies are shown using circles. 



explained by 'small accretion events' which contribute 
blue members to the clusters without much increase of 
other parameters such as mass or X-ray luminosity. Such 
events could be a source of scatter in the galaxy popu- 
lations of clusters selected by any selection method. To 
analyze correlatio n with other cluster parameters |De| 



Propris et al. ( 2004 1 used a sample of clusters from 2dF 
Galaxy RedskiirSurvey (2dFGRS) at redshift < 0.11, 
finding a large variation (/;, ^ 0.1-0.5 for Al* 4- 1.5 at 
'"200) from cluster to cluster. 

The apparent contradiction between X-ray and opti- 
cally selected samples and the sensitivity to environmen- 
tal effects, raises questions about how much of the ob- 
served fb is due to a selection method, how much it is 
due to the intrinsic scatter, and if these two effects can 
conspire to produce an apparent trend where no trend 
exists. 

What is needed is a sample of galaxy clusters which 
possess two main characteristics: (1) the selection of clus- 
ters is made in a way that is independent of the quantity 
whose evolution is being studied to avoid pos sible bias 
( [Newberry et al.||1988[ [Andreon fc Ettori||1999[ ), and (2) 
the sample must contain the same class of clusters (i.e. 
same mass range) at different redshift to h elp in sepa- 



rating mass trends from redshift evolution (Andreon & 
| Ettori|1999 ) . A sample of SZE selected clusters of galax- 
ies fulfills tnese requirements. The selection of the SZE 
clusters is closely related to mass, and that mass selec- 
tion is approximately independent of redshift, allowing 
a comparison of the same type of clusters at different 
epochs. 

Historically fb has been measured in different ways. 
Initially the average color of the E/SO galaxies, within 
a radius of R30 from the cluster center that is the ra- 
dius that contains 30% of all galaxies that belong to the 
cluster, and a conc entration index, were use to define fb 
(sec Butchcf & Ocmlcr 1984; Rakos fc Schombcrt' |1995' 
|Margonincr fc dc Carvalho.2000^ I'airley et al.^2002| ). An- 
other approach is using the red sequence from the color- 
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masnitude d i agram of the clusters and r2oo (Popesso 



et al. 



2007b 
ot both met 
agram but R30 



Barkhouse et al. 2007) or a combination 



rods , that is using the color-magnitude di 



20021 



(Kodama & Bower 2001 Fairley et al 



Here we follow the approach of using the red sequence 
to define the red and blue populations, and r2oo to de- 
fine the radial extent. This ensures we are using the same 
portion of the cluster virial region, independent of red- 
shift, and that we are exploring populations with colors 
defined with respect to a passively evolving SSP model. 

The galaxies used for the ft measurement are inside 
the r2oo cluster radius and are fainter than the BCG and 
brighter than 0.36L*j^qjj^j^. They are classified as red 
if they are located within ±3 times the average disper- 
sion of the G aussian fit to the color-m agnitude relation 
(± 0.22 mag; Lopez-Cruz et al. 2004), and blue if they 
are more thaii 0.22 mags bluer than the Red Sequence. 
We choose a limit of 0.36L*j^qjj^j^ to allow a meaning- 
ful comparison among three of our four clusters, as it 
is the deepest magnitude that we can detect with good 
completeness for the three of them. For the fourth clus- 
ter, SPT-CL J0546-5345, we currently do not have deep 
enough photometry for this analysis. 

The color magnitude diagram used for the clusters de- 
pends on the red sequence identification: g-i/i for SPT- 
CL J0516-5430, r-i/i for SPT-CL J0509-5342, and i-z/z 
for SPT-CL J0528-5300 and SPT-CL J0546-5345 (see fig. 
11 ). The blue fraction is defined as the statistically back- 
ground corrected number of blue galaxies rif, divided by 
the total number of statistically background corrected 
galaxies nt. The blue fraction and its gaussian propa- 
gated uncertainty are: 



fb 



rib 



E 



9/ 



(1) 



Where rib and are the blue and red statistically back- 
ground subtracted number of galaxies: 



The uncertainties are expressed as 



ibkg) 



assuming aj^. Poissonian. The last term is calculated 
directly by measuring the RMS of the Gaussian distribu- 
tion observed on histograms constructed from the blue 
and red (or total) number of galaxies background cor- 
rected in a circle of radius r2oo on n random position out- 
side the cluster radius (ha,ckgTomid(r2oo)-i>ackground) in 
order to account for background variations on the ob- 
served 36'x36' patch of the sky. 

A special mention for SPT-CL 0509-5342 is required. 
In the center of the cluster are three bright stars leav- 
ing only a few visible galaxies; we have corrected this 
effect by accounting for the area masked around these 
stars. Nevertheless, the statistical background subtrac- 
tion leads to negative blue galaxy counts in the ~ 0.6r2oo 
inner part of the cluster area. 

The blue fraction of three of the four clusters, at 
redshifts 0.295, 0.463 and 0.763, are shown in Fig. ^ 
The measurements suggest an increase with redshift, as 
shown for optically selected clusters, although the result 
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Fig. 12. — Blue f ract ion ve rsus redshift using the po pula- 
tions shown in Figure [Tl] The [Butcher fc OemIer| | |1984| l rela- 
tion is shown (dashed line) . The blue traction is calculated using 
statistical background correction within r2oo and to ^ depth of 
'^(>L'^jODEL each system. 

could be consistent with a constant blue fraction over the 
range of redshift that we explored with our limited sam- 
ple. Future optical follow up of SPT-SZE selected clus- 
ters using larger aperture telescopes on the high redshift 
end will be necessary to understand the Butcher-Oemler 
effect in this cluster mass range within the SZE selected 
sample. 

5. CONCLUSIONS 

We present the results of a careful examination of the 
multiband optical properties of the galaxy populations in 
the first four SZE selected galaxy clusters. This analy- 
sis builds upon the selection by the South Pole Telescope 
survey, deep multiband optical data from the Blanco Cos- 
mology Survey, Chandra and XMM mass estimates and 
published spectroscopic redshifts. 

The radial distributions of galaxies in the four systems 
are consistent with NEW profiles with low concentra- 
tion in the 2.2-3.6 range, although the constraints in 
our highest redshift clusters are weak due to the imaging 
depth. One system shows a clear secondary peak, which 
is evidence of multiple galaxy components. The observed 
galaxy concentrations in these SPT systems are consis- 
tent with X-ray and optical selected cluster samples as 
well as simulations. 

We showed that the characteristic luminosities in 
bands griz are consistent with passively evolving pop- 
ulations emerging from a single burst at redshift z = 3. 
This is observed by direct comparison of the griz m* 
measurements with the evolution of the red sequence ex- 
pected from the SSP model. 

The slope of the luminosity function, a, in all four 
bands showed an average of —1.2 consistent with previ- 
ous studies and roughly independent of redshift, although 
in the high redshift systems the constraints are weaker 
and the $ — a contours are much more extended (see 
Fig. [9]) due to the depth of the data. 

Halo occupation numbers (to m* + 2) for these sys- 
tems appear to be consistent with the relation measured 
in X-ray s elected clusters. As shown previously (Lin 
et al.|2004 1 , this well behaved and simple galaxy popula- 
tions IS unfortunately not easy to use as a mass indicator 
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with optical data alone, because the HON varies with the 
adopted virial radius of the cluster. 

The blue fractions observed in these systems are con- 
sistent with those seen in clusters selected using other 
means. Although the meassured fb suggest a redshift 
evolution (as optical studies show), it is within the errors 
also consistent with a constant ft,. It is clear that defini- 
tive conclusions should be drawn with a larger number 
of clusters for more robust statistics. A larger sample 
and deeper multiband data on the high redshift end is 
needed. 

The SPT selection provides a powerful means of choos- 
ing similar mass systems over a broad range of redshift, 
making the future larger cluster sample particularly in- 
teresting for this study. 

In summary, our systematic analysis of the galaxy pop- 
ulations in the first SZE selected galaxy clusters spanning 
the redshift range 0.3 < z < 1.1 provides no clear evi- 
dence that the galaxy populations in these SPT selected 
clusters differ from populations studied in other X-ray 
and optically selected samples. An extension of our anal- 
ysis to the full SPT sample will enable a more precise 
test of the effects of selection. In addition, comparison 
of the observed properties of the SPT cluster galaxy pop- 
ulations and their evolution to numerical simulations of 
galaxy formation should allow for clean tests of the range 
of physical processes that are responsible in determining 



the formation and evolution of cluster galaxies. 
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